rm(list=ls())

# Figure A.1

exp_claassen_m5 <- readRDS("claassen_mod5_output.rds")

pdf("fig_a1.pdf", width=10, height=4)
rstan::traceplot(exp_claassen_m5, ncol=5, nrow=2, alpha=0.8, size=0.3, 
                 pars=c("mu_lambda","sigma_lambda","sigma_theta","sigma_delta","phi","lambda[15]",
                        "delta[107]","theta[21,66]","theta[16,49]","theta[16,101]"))
dev.off()
